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In this paper, we show that there is a large class of fermionic systems for which 
it is possible to find, for any dimension, a finite closed set of eigenoperators and 
eigenvalues of the Hamiltonian. Then, the hierarchy of the equations of motion 
closes and analytical expressions for the Green's functions are obtained in terms 
of a finite number of parameters, to be self-consistently determined. Several 
examples are given. In particular, for these examples it is shown that in the 
one-dimensional case it is possible to derive by means of algebraic constraints 
a set of equations which allow us to determine the self-consistent parameters 
and to obtain a complete exact solution. 
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One of the most intriguing problem in Condensed Matter Physics is the study 
of highly interacting fermionic systems [[1]. The standard methods based on pertur- 
bation theories do not work well and the attention is put on developing alternative 
approximate formulations. In this context, it is interesting to cultivate the study of 
some solvable models, which are themselves of physical interest and which may fur- 
nish some indication on the solution of more complex models. In a recent paper [[2], 
we have shown that there is a large class of fermionic systems for which it is possible 
to find a complete set of eigenoperators and eigenvalues of the Hamiltonian. Then, 
the hierarchy of the equations of motion closes and analytical exact expressions for 
the Green's functions (GF) can be obtained. 

We consider a system of q species of particles, satisfying Fermi statistics and 
localized on the sites of a Bravais lattice. We suppose that the mass of the particles is 
very large and/or the interaction is so strong that the kinetic energy is negligible and 
the particles are frozen on the lattice sites. Let Ca{i) and cl{i) be the annihilation and 
creation operators of the particles of species a in the Heisenberg picture: i = (i, t), 
where i stands for the lattice vector Ri. These operators satisfy canonical anti- 
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commutation relations 

{Ca(i,t),4(j,t)} = 5afe5ij /^v, 

K(i,t),c,(j,t)} = {ct(i,t),4(j,t)} = ^ ^ 

A quite general Hamiltonian for these systems reads as 

ia ij ab ijl o-bc 

(2) 

where Vaii) represents an external field acting on the particle a; Uaii) = cl{i)ca{i) 
is the particle density operator of the species a; Vat, Vahc , ■ ■ ■ are the two-body, 
three-body, . . . potentials. In general, the potentials will be translationally invariant 
and will depend only on the differences between the coordinates of the particles. 

The class of systems described by the Hamiltonian ([2]) is very large and of direct 
physical interest. By appropriate choice of the potentials it may describe multiband 
electronic systems dominated by charge correlations and by magnetic interactions 
in the s^-channel, where is the third component of the spin. Furthermore, let 
^(0 = Z]acl(Oca(^) be the total particle density and consider the transformation 
n{i) = 2 + 'S'(^)- Under this transformation the Hamiltonian ([2]) is mapped to spin- 
q/2 Ising-like models. 

Now, if the interaction potentials have a finite range, the model Hamiltonian is 
always solvable. The proof of this statement is the following. Firstly, it is immediate 
to see that the density operator na{i) does not depend on time 

i^^naii) = [na{t),H] = (3) 

Then, in order to use the equations of motion and GF formalism we must start from 
the Heisenberg equation for the fermionic field. It is immediate to see that Ca{i) 
satisfies the equation of motion 

d 1 

i oTCa(«) = Va{i)Ca{i) + ^ j)^6(j)Ca(«) + ^YYl ^abc{i,i, l)n;,(j>c(/)Ca(i) + • • • 

jfe jl be 

(4) 

The dynamics has generated other field operators of higher complexity. By taking 
time derivatives of increasing order, more and more complex operators are generated. 
These operators might be named composite operators, as they are all expressed in 
terms of the original fields Ca{i)- Because of the finite range of the interaction, the 
higher-order field operators generated by the dynamics extend over a finite number 
of sites. Then, because the particle operator satisfies the algebra [na{i)Y = ^aii) 
(for p > 1), the number of composite operators is finite and a complete set of 
eigenoperators of the Hamiltonian can be found. If n is the number of independent 
composite operators, we can construct a n-multiplet operator 

/ ^iW \ 

Mi) 



(5) 
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which satisfies the Heisenberg equation 

i^^m = [m.H] = em (6) 

where the n x n matrix e will be denominated as the energy matrix. Note that we 
are using a vectorial notation: the field ipmii) is itself a multiplet of rank q. Once the 
composite operator and the energy matrix have been determined, an exact solution 
of the Hamiltonian can be formally obtained. Let us define the retarded Green's 
function 

Gi^,J) = (Rimi^^j)]) = e{U - ({^(z), (7) 

where (■ ■ ■) denotes the quantum-statistical average over the grand canonical en- 
semble. By means of the Heisenberg equation ([S]), we obtain in momentum space 
the equation 

[cj-e]G(k,cj) = /(k) (8) 
where /(k) is the Fourier transform of the normalization matrix, defined as 

/(i,j) = ({^(i,t),^t(j^i)|^ (9) 

The solution of Eq. ([8]) is 

^-^uj - Em + lb 

where Em are the eigenvalues of the energy matrix e. The spectral density matrices 
(j{m) ^ jj.^ g^j,g calculated by means of the formula [ [3] 

ai'^)(k) = (]«™$:[fi^,]-%^(k) (11) 

7 

where Vt is the matrix whose columns are the eigenvectors of the matrix e. The 
correlation function C{i,j) = can be immediately calculated from ( |TOl) 

and one obtains 



C(k, uj)=7r YTmli - ^™)T™a(™) (k) 



(12) 



with Tm = 1 + tanh (/3i?m/2), /3 is the inverse temperature. By similar technique 
we can easily calculate [H] multi-point correlation functions as C{i,j; li, h-, - ■ ■ h) = 

Equations (|T0|) and ([T2l) are an exact solution of the model Hamiltonian ([2]). 
However, the knowledge of the GF is not fully achieved yet. The algebra of the field 
is not canonical: as a consequence, the normalization matrix /(k) in the equation 
of motion ([8]) contains some unknown static correlation functions. Generally, these 
correlators are expectation values of operators not belonging to the chosen basis 
and should be self-consistently calculated. According to the scheme of calculations 
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proposed by the composite operator method [ [3l [5], one way of calculating the 
unknown correlators is by specifying the representation where the GF are realized. 
The knowledge of the Hamiltonian and of the operatorial algebra is not sufficient 
to completely determine the GF. The GF refer to a specific representation (i.e., 
to a specific choice of the Hilbert space) and this information must be supplied to 
the equations of motion that alone are not sufficient to completely determine the 
GF. The procedure is the following. From the algebra it is possible to derive several 
relations among the operators. We will call algebra constraints (AC) all possible 
relations among the operators dictated by the algebra. This set of relations valid at 
microscopic level must be satisfied also at macroscopic level, when expectation values 
are considered. Use of these considerations leads to some self-consistent equations 
which will be used to fix the unknown correlators appearing in the normalization 
matrix. An immediate set of rules is given by the equation 

{mi^\^)) = ^j:^Jd^C{k,u;) (13) 

k 

where the l.h.s. is fixed by the AC and the boundary conditions compatible with 
the phase under investigation, while in the r.h.s. the correlation function C(k, tu) is 
computed by means of the expression f|T2l) . Use of f|T3|) will lead to a set of exact 
relations among the correlation functions. Unfortunately, for the considered class 
of systems, the number of equations is always less than the number of correlation 
functions and extra equations are needed. To this purpose we have recently found [ 
[21 H] a set of AC capable to give other relations among the correlation functions. Let 
us suppose that there exist some operators, O, which project out of the Hamiltonian 
a reduced part 

OH = OHo (14) 

When Ho and Hj = H — Hq commute, an important rule which descends from the 
AC (1141) is that the quantum statistical averages over the complete Hamiltonian H 
must coincide with the average over the reduced Hamiltonian Hq 

Tr{Oe-^"} = Tr{Oe-^^°} (15) 

For a variety of one-dimensional models, we have found that use of the condition 
ffT^ allows us to close the set of equations for the correlation functions and to obtain 
an exact solution of the model. For dimensions higher than one these conditions are 
not sufficient and more equations are necessary. 

We now apply the above procedure to some specific models, obtained by special 
choice of the interaction potentials. We shall consider a d-dimensional hypercubic 
Bravais lattice. 

First example 

We consider a g-state model described by the Hamiltonian 

= -/i^n(z) + c/V^n(i)n"(i) (16) 
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where fi is the chemical potential, n{i) = J2a^a{i) is the total particle density. 
Hereafter, for a generic operator $(i) we us the notation <l>"(i) = X^j ttij'^'O, t), where 
Oij is the projector on the first-nearest neighbor sites. This model might describe 
a system of g-electron bands, interacting through an intersite Coulomb potential 
of strength V. Alternatively, by means of the transformation n{i) = f + S{i) the 
Hamiltonian flTB]) takes the form 

H = -hY^S{i) - + Eo (17) 

i i 

where J = —dV, h = fi — qdV, Eq = |(— yU + ^dV)N; N is the number of sites. 
Hamiltonian f[T7|) is just the spin-g/2 Ising model with nearest neighbor interactions 
in presence of an external magnetic field h. The following recursion rule can be 
established [3] for the operator [?2"(i)]P 

2qd 

K(.)f = E A(^)K(0]- (18) 

m=l 

where the coefficients are rational numbers, satisfying the relation J2m=i ^m'* ~ 
1, that can be easily determined by the algebra and by the structure of the lattice [ 
H]. Owing to this rule, the set of eigenoperators of the Hamiltonian is given by 



( c(0 \ 

c{i)n°'{i) 

V c(i)[n"(i)]2<?'^ / 



{i),H] = em (19) 



where the energy matrix e can be calculated by means of the equation of motion 
and the recursion rule (ITS]) . The eigenvalues E^ of e are given by 

E^ = -^ + {m-l)V (m = l,2,---2gd+ 1) (20) 

The retarded Green's functions and the correlation functions can be exactly calcu- 
lated by applying the scheme of calculations illustrated above. Then, by using the 
AC ( TTH) it is possible to derive the self-consistent equations 



2qd+l 



^(.-1) _ ^(.-1) ^ ^ ^ y^^M (A; = 1, . . . , 2gd + 1) (21) 

^ m=l 

where Tm = 1 + tanh (/3Em/2) and a^'^^ are the spectral density matrices, that can 
be calculated by means of ( ITTl) and are expressed in terms of the quantities 
Equations fl2T]) give a set of exact relations, valid for any value of q and d, among 
the correlation functions 

«:(p) = (K(On 

X(P) = {n{z)[n'^{z)r) (p-0,l,---2gd) (22) 
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Unfortunately, the number of equations is not sufficient and we need other condi- 
tions. We have studied the cases of g = 1 ,2, 3 in Refs. [SI El [7], respectively. In 
all these cases for one-dimensional systems, it is possible to find by means of the 
algebraic constraints f|T5|) the necessary equations to close the set. These extra con- 
ditions are obtained by exploiting the following algebraic conditions. Let us divide 
the Hamiltonian as 

H = Ho + 2Vn{'i?)n''{i) (23) 

where, due to translational invariance, i is a generic site of the infinite chain. Then, 
we have 

q = 1 ct(«)e-^^ = cH2)e-^^° (24) 
^^{i)e-(^^ = e^(z)e-^^» 
^ " ^ r]^{i)e-^^ = r]^{i){l + E /^K(i)]'"}e-^^° ^^^^ 

m=l 

q = 3 vK^)e~^'' = vK^){l + /«^K(0]"^}e-^^° (26) 
For g = 2 the definitions are: 



CKi)e-^" = C^(0{1 + E (2/™ + ^7™)K(z)]-}e-^^o 

m=l 



^a(i) = [1 - n{{}]ca{{} .27) 
Va{i) = n{i)ca{i) 

The operators C,a{i) and ria{i) induce the transitions -vv- 1, 1 -v^ 2, respectively. For 
g = 3 the definitions are: 

U^) = [^-n{^) + D{^)]ca{i} 

r)a{i) = [n{^ - 2D{^]ca{^ (28) 
Ui) = D(t)c,(t) 

where D{i) is the double occupancy operator, defined as 

D{{) = ni{i)n2{i) + ni{i)n^{i) + n2{i)n^{i) (29) 

The projection operators i^a, and C^a induce the transitions among states with 
different particle numbers: <S=v- 1, 1 -vv- 2, 2 <;=^ 3, respectively. The quantities fm 
and Qm are known functions of 13V . 
Second example 

Let us consider two species of particles, say a and 6, and consider the Hamiltonian 
H = -fiJ2 ^(^) + UJ2^('^ + n(2)n"(i) (30) 

i i i 

where n{i) = na{i) + nh{i) and D{i) = na{i)ni,{i) are the total particle density and 
double occupancy operators, respectively. This Hamiltonian is just the extended 
Hubbard model in the ionic limit, where U and V are the on-site and inter-site 
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Coulomb interaction, respectively. The two species of particles, a and b, are in this 
case electrons with spin up and down, respectively. By means of the transformation 
n{i) = 1 + S{i) , (|30|1 can be cast in the form 



(31) 



where J = -dV, h = fi- 2dV - ^U, A = ^U, Eq = + dV)N. Hamiltonian §^ 
is just the Ising spin-1 model with nearest-neighbor interactions in the presence of 
a crystal field A and an external magnetic field [[HI [9l [lOl [TT]. We now define the 
composite operators 



) \ 



Ml 



r]{i 



) \ 



(32) 



V r7(i)[n"(i)]4^ / 



where ^(i) and ri{i) are the Hubbard operators, defined in (p71) . By means of (fT8|) . 
these fields are eigenoperators of the Hamiltonian (!30|) 



(33) 



where e^^-* and e^^-* are the energy matrices, of rank (Ad + 1) x (Ad + 1) , which can 
be calculated by means of the equations of motion and the recursion rule (fTSjl . The 
eigenvalues and of the energy matrices are given by 

Ei^ = -f, + (m-l)V 
E^2^ = -^ + U +(m~l)V 



[m 



1,2,---Aqd + 1] 



(34) 



The retarded and the correlation functions can be exactly calculated by applying the 
scheme of calculations illustrated above. Then, by using the AC fll3p it is possible 
to derive the self-consistent equations 



2 



1 Ad+i 



(k 



■Ad+l) (35) 



m=l 



where = 1 + tanh (^f3E^^ with a = ^,77. o"*^"'™^ are the spectral density 
matrices, that can be calculated by means of ffTTj) and are expressed in terms of 
the quantities k^^^ and A'^^^ Equations fl35l) give a set of exact relations, vahd for 
any dimension, among the correlation functions fl2^ . Unfortunately, the number of 
equations is not sufficient and we need other conditions. Our study [ E] has shown 
that for the one-dimensional case, it is possible to find by means of the algebraic 
constraints (fT5|) the necessary equations to close the set (1351) . These extra conditions 
are obtained by exploiting the following algebraic conditions 



l3Ho 



^^(t)e-^H = ^t(z)e- 
D(t)e-f'' = D(t){l+ E (2/^ + (?^)[n-(z)]-}e-'^^° 

m=l 



(36) 
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where = H — 2Vn{i)n"{i). 
Third example 

The model fl30l) can be generalized by considering 3- and 4-body potentials. In 
particular, by considering one-dimensional systems, let us take 

i i i i i 

(37) 

This Hamiltonian can be mapped to the following model 

H = -hJ2 S{t) + A 5] S\t) - - KY^S^^jS^'^ii) + Eo (38) 



where J 



(— /X + V)N. This Hamiltonian is just the one-dimensional Blume- Emery-Griffiths 
model [[T2]- We now define the composite operators 



( m \ 

e(^)K(^)]2 

e(^)[n-(2)]3 



( \ 

r]{i)rf{i) 
?7(i)[n"(i)]2 
r7(i)[n"(i)]3 
?7(i)[n°(i)]^ 
r/(i)L)°(i) 



(39) 



where ^{i) and rj{i) are the Hubbard operators, [cfr. (1^ ]. These fields are eigenop- 
erators of the Hamiltonian 0371) . The corresponding eigenvalues are 



^(0 



-^ + v 

-fi + 2V 
-fi-U + 2V 
-fi-^U + 2V 
-fi-U + AV 
[-fx-^U + 3V I 



I 



hu + v 



— jj^ — ^ — 

-ti--f-U + 2V 
-fi--f + 2V 
-\U + 2V 
-U + AV 
-fi--f~U + 3V J 



(40) 



The retarded and the correlation functions can be exactly calculated by applying the 
scheme of calculations illustrated above. Then, by using the AC f|T3|) it is possible 
to derive [[13] the self-consistent equations 



(A; = l,---5) 

(k = 6, 7) 



(41) 



where the definitions are the same as in example 2, and the new correlations functions 
and are defined as 



S(p) = ([D°(i)]P) 
^(p) = (n(i)[D"(i)]P) 



(P=l,2) 



(42) 
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the spectral density matrices and are expressed in terms of the quantities 
^(p)^ \ip)^ §ip)^ gip)_ Equations fHT]) give a set of exact relations among the correlation 
functions fl22|) and fj42l) . By exploiting the algebraic condition C,^{i)e~^^ = C,^{i)e~^^° 



it is possible to derive 

'-'1,4 
'-'1,7 

where 



the self-consistent equations 



Cjfi [^Xi + X2 + ^Xf] 

^1,1 [i^l + 1^2 + 1-^1-^2 + 

cjfi [IX2 + 



(43) 



'-'1,2 

MO 
'-'1,1 



X, 



CI 



(0 



MO 
'-^1,1 



(44) 



C[^l is the equal time correlation function C[^l = expressed in 



(45) 



terms of the parameters S^^\ 9^^^ by means of the relation 



^ m=l 



k 



Fourth example 

Let us consider the case of two particles, characterized by spin cr =1, | 
and choice the potentials as 

K-(i) = — /U — ah 

K,<,.(i,j) = 2^6/(2(5,,,, -l)«ij 

h is the intensity of the external magnetic field. The Hamiltonian is 



(+,-), 



(46) 



(47) 



where n^{i) = n-i{i) —ni{i) is the third component of the spin density operator. Let 
us restrict the analysis to one-dimensional systems. The following recursion rule can 
be established for the operator [^3 (i)]^ 



(4^ 



m=l 



where the coefficients A^p^ are rational numbers, satisfying the relation A^p^ = 

1, that can be easily determined by the algebra and the structure of the lattice. 
Owing to this rule, the set of eigenoperators of the Hamiltonian (H7!) is given by 




c„{i)n%{i) 
\ c.(z)[nf(z)]^ / 



(49) 
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(50) 



The energy matrix e^'^^ is a 5 x 5 matrix with the following expression 



f —fi — ah 2aV 

— yU — ah 





V 









2aV 
-fi — ah 





The eigenvalues of the matrix e^"") are 



—fi — ah 
-l2-ah-2V 
—fi — ah — V 
—fi — ah + V 
\ -H - ah + 2V J 






2aV 
-fi — ah 
\aV 







2aV 
-fi — ah J 



(51) 



(52) 



The retarded Green's functions and the correlation functions can be exactly calcu- 
lated by applying the scheme of calculations illustrated above. Then, by using the 
AC ( IT3l) it is possible to derive the self-consistent equations 



1 y^5 



is 



=1 [^m ■* 



a 



l,k 



I 7TJ. 



a 



(T,m) _ 



l,k 



l,k J 



(53) 



where Tj^^ = 1 + tanh (^pE^^^ /2) and a^'"'") are the spectral density matrices, that 

can be calculated by means of (JTTll and are expressed in terms of the quantities 
Equations fl53l) give a set of exact relations among the correlation functions 



Xip) 
0{p) 



([n^(z)n 

(n(z)[n-(0]P) 

{n,{i)[nt{i)Y) 



(p = 0,l,---4) 



(54) 



Also for this model we can use the AC f|T5l) in order to derive extra equations capable 
to close the set fl53l) for the correlation functions. Details will be presented elsewhere. 

Summarizing, we have shown that a large class of systems of localized particles, 
satisfying Fermi statistics and subject to finite-range interactions, is always solvable, 
in the sense that a complete finite set of eigenoperators and eigenvalues of the 
Hamiltonian can be found. This knowledge allows us to derive analytical expressions 
for the Green's functions and for the correlation functions and a set of exact relations 
among the correlation functions can be derived. As an illustration we have considered 
several examples. In all the studied models we have shown that for the case of one 
dimension it is possible to use algebraic constraints which permit to close the set of 
self-consistent equations and to obtain exact solutions. For higher dimensions more 
self-consistent equations are needed. This problem is now under investigation. 
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